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We revisit the problem of phantom behaviour of effective dark energy in scalar- 
tensor gravity. The main focus is on the properties of the functions defining the 
model. We find that models with the present phantom behavior can be made con- 
sistent with all constraints, but one of these functions must have rather contrived 
shape, and the initial data must be strongly fine-tuned. Also, the phantom stage 
must have begun fairly recently, at z < 1. All this disfavors the effective phantom 
behaviour in the scalar-tensor gravity. 



I. INTRODUCTION 

Most models of dark energy in the present Universe predict that its effective equation 
of state satisfies the null energy condition (NEC) w e ff = Pde/pde > — 1, where pde 
and pde are the effective dark energy density and pressure, respectively. However, the 
observations do not rule out that dark energy is phantom, i.e., it violates NEC. As 
an example, the 7-year WMAP+BAO+SN data pQ give the following bound on the 
equation of state with time-dependent w e ff at z = : 

Weff = -0.93 ± 0.13 (68%C7L) , 

which is not entirely inconsistent with iu e //,o < — 1- Even though phantom dark energy 
can be accomodated within General Relativity [2H5], it is legitimate to ask whether ef- 
fective phantom behavior can be obtained in modified gravity theories, such as f(R) 
or scalar-tensor gravity [6HTT]. This question can be addressed, in particular, within 
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a popular approach [6j [121 H3] employing the reconstruction of model parameters from 
the redshift expansion of observables, combined with the experimental constraints on 
non-GR gravity. The main conclusion is that with current observational data, the recon- 
structed DE behaves more or less like the cosmological constant, but it is still possible 
to have phantom stage today. 

In this paper we follow a somewhat different route and ask what sort of the scalar- 
tensor Lagrangian can lead to effective phantom DE at the present epoch without vio- 
lating constraints from Solar system and local gravity experiments, time- (in) dependence 
of the gravity constant, etc. We also ask whether fine-tuning of initial data is necessary 
and how long the duration of the phantom stage can be in the past. 

Our results are somewhat disappoining. We find that models with the present phan- 
tom DE can be made consistent with all constraints, but one of the functions entering 
the scalar-tensor Lagrangian must have rather specific shape, and the initial data must 
be strongly fine-tuned. Also, the phantom stage must have begun fairly recently, at 
z < 1. Before that the scalar field was undistinguishable from quintessence. All this dis- 
favors the effective phantom behaviour in the scalar-tensor gravity. In fact, as we point 
out towards the end of this paper, some of the unpleasant properties we discuss must 
be present in scalar-tensor models for effective dark energy irrespectively of whether it 
is phantom or not. 

This paper is organized as follows. In Section II we present the equations governing 
the homogeneous cosmological evolution in the scalar-tensor gravity. We recall in Sec- 
tion III the experimental constraints on the non-GR gravity, that place bounds on the 
parameters of the theory. In Section IV we define the expansion coefficients of the func- 
tions entering the Lagrangian and reformulate the bounds of Section III in terms of these 
coefficients. In Section V we put together all consistency requirements for the effective 
phantom dark energy today and arrive at qualitative understanding of the properties of 
the functions defining the theory. Also, the maximum redshift at which the phantom 
phase could begin is estimated. Section VI contains a numerical example. We conclude 
in Section VII. 
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II. HOMOGENEOUS AND ISOTROPIC EVOLUTION 

By definition, the effective dark energy density and pressure p e ff and p e jf = w e ffp e ff 
are the quantities entering the GR-looking evolution equations for the homogeneous and 
isotropic Universe, 

3H 2 = p + p eff (1) 
-2H = p + p eff (l + w eff ) (2) 

where p and p = are matter energy density and pressure, and we set 87tGn = 1- Using 
and @, one writes 

where Q m = p/3H 2 . We are going to make use of this relation in the context of the 
scalar-tensor gravity. The action of this theory is 



S = \ j d 4 x^ (F($)R - Z^g^d^d^ - 2C/($)) + S m (^, 9^), 



(4) 



(mostly positive signature), where the action for the usual matter S m does not depend 
on $. One can always redefine the field to have a convenient form of either or 
Z($). We will use the general form of F(§) and set 

Z($) = 1. 

From the action ^ one obtains the gravitational equations, 

(r^-^r) =T^ + d^d^- l -g liy {d^) 2 + V l N v F^)-g^UF^)-g ilv U^) 

(5) 

and equation of motion for the field $, 

ldF dU 
2 d<3> d$ 

Let us specify to the homogeneous, isotropic and spatially flat Universe with metric 



ds 2 = -dt 2 + a 2 (t)cfeW. 
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Since matter does not interact with $, the scale factor a has the same meaning as in 
GR. Using Eqs. ^ and ^ one gets the following set of equations: 

3FH 2 = p + ^<i> 2 - 3HF + U (7) 

-2FH = p + <5> 2 + F - HF (8) 

, ■ 0n dF dU , , 

$ + 3#$ = 3 (H + 2H 2 — -— . 9 

The equation for the matter density has the usual form, 

p + 3Hp = 0. 

It is convenient for futher analysis to switch from the evolution in time to the evolution 
in redshift. This can be done by using the relation 

In this way one obtains from ([7]), ^ and ^ the evolution equations in terms of redshift: 

3FH 2 = p + H 2 (l + z) 2 — + 3H 2 (l + z)F' + U (10) 

2FHH' = p + H 2 (l + z)V 2 + H 2 (l + z) 2 F"+ 

+ [(1 + z) 2 HH' + 2H 2 {1 + z)]F' (11) 

# 2 (1 + z) V + [M'(l + ^) 2 - 2H 2 (l + z)]& = 3[2H 2 - HH'(1 + z)]^ - ^, (12) 
where prime denotes d/dz. 



III. CONSTRAINTS 

The properties of functions defining the theory are strongly constrained by local and 
Solar system experiments. This is a major problem for the effective NEC-violating 
behavior in the scalar-tensor gravity. One important parameter is the Brans-Dicke 
"constant" Wbd{&)- It is straightforward to obtain the expression for Wbd m our 
parametrization by redifining the scalar field. One finds 



Wbd = — (13) 

\d<5>) 
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The lower bound on the present value of Wbd is obtained from the Kassini experiment 
(HHTE]. It reads (the subscript denotes the quantities at the present epoch) 

W BDfi > 4 • 10 4 . (14) 

A bound of another sort follows from the experiments on the time dependence of the 
gravity constant. In our case the local gravity constant G\ oc is given by |Bj 

_ 1 ( 2F + A{dF/d^f \ _ l 2W BD + 4 
loc ~ F \2F + Z{dF/d<5>Y ) ~ F2W BD +3' 1 J 

The experimental constraint on the time evolution of Gi oc can be found in Refs. [TTJUE]. 
For z — it reads 



HG 



loc 



<0.5-l(T 2 . (16) 
o 



We also know that the gravitational constant relevant for cosmology should not change 
significantly since Big Bang Nucleosynthesis [19J, 



^<0.1. (17) 

^cosm 

It is the latter constraint that plays a significant role in our analysis, see Section V. 



IV. EXPANDING IN REDSHIFT AND $ 

A convenient way to analyze the evolution at small redshifts is to expand all functions 
in the Taylor series in redshift z. On the other hand, we are mainly interested in the 
dependence on $, so we will use the mixed expansion. At z = 0, without loss of generality 
we choose 

$o = 

and by definition of the Newton gravity constant we have 



F = 1. 
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Here is our definition of the expansion coefficients: 



F(z) = 1 + F x z + \f 2 z 2 + \f z z\.. 

2 o 

u{z)/ml = n Ufi + u lZ 

H 2 (z)/H 2 = l + h lZ +h 2 z 2 + ... 
V{z) = %z + U>lz 2 

p(z)/3h 2 = n mfi (i + zf. 



(18) 
(19) 
(20) 

(21) 
(22) 



Without loss of generality we take $' > 0. From Eq. (fO) it is straightforward to obtain 



the relation between the derivatives at the present time, 



$' 2 = 6(1 -n u>0 -n m , -Fx). 



(23) 



We will also need h\ and h 2 to obtain the expression for u> e //. From Eqs. (10), (11) we 
get 



hi 

h 2 



1 



1 - 



K (6 - 3fi m , - 60^,0 - 4Fi + F 2 



(24) 



(1-f) 

+ -^l ^—"17,0 



2 [F x ( |f x - 3^F 2 



2 3F2 o 

— — Sir/n — 



2 2 
3F 2 



fi m>0 + 2F 2 + — - 5fi(7 ; o - mi - 4f2 mj0 + 5]. 
4 6 



(25) 



Our main purpose is to understand the behavior of F and U as functions of the scalar 
field, so we expand them in 



F = 1 + /i$ + ^/ 2 $ 2 + ^/ 3 $ 3 , 

1 9 

£7 = 1/0 + "I* + 2 M 2$ • 



(26) 
(27) 



The relationship between the expansion coefficients entering (18) and (26) is 



Fx = fi&o 

F 2 = fl$o + f2% 

F 3 = /!< + 2f 2 $»& + / 3 $ /3 . 



J2 



(28) 
(29) 
(30) 
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We now recall the constraint on Wbdo, Eq. (14), and make use of Eq. (13). With 



our normalization F — 1, we get very strong upper bound on /j.: 



\h\ <0.5-l(T 2 . 



(31) 



This means that the field $ is presently near the extremum of the funcion F{<&). Such 
a conclusion appears inevitable in modified gravity, see, e.g., Ref. [TJ. 



It follows from Eq. (23) that $' 2 < 1, so Eq. (28) implies that F\ is also small 



l^il < io- 2 . 



(32) 



This suggests that we can neglect terms with the first derivative of F in the analysis of 
the present epoch. We note in passing that a small value of F'(z = 0) could have been 
anticipated, since the GR tests are very precise, and only a slight deviation from GR 
can be tolerated today. 



Using (23) and neglecting the term with iq, we obtain for the present value of the 



field derivative with respect to redshift 



$^ = 6(1-^,0 -ft 



m,0J 



(33) 



This simple relation will be instrumental in what follows. 



V. EFFECTIVE PHANTOM BEHAVIOR 



Now we use the expression ^ for ti> e //,o to find out which parametrs can be re- 



sponsible for the effective phantom behavior. Making use of Eqs. (28), (29) and (32) we 
get 

, ^ + 6(1-^,0-^,0) + /!^ 



3(1 -tt 



m,0) 



By extracting the second derivative of the field from Eq. (12), we find 

hi 



< - U x . 



So, there are essentially two parameters that could yield 1 + w e //,o < 0, namely, f 2 and 
U x . We begin our discussion with the latter. 
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The possible contribution of the potential to the phantom effective equation of state 



comes from the third term in (34) (the second term in the numerator is positive in virtue 



of Eq. (33)) and is given by 

At,(l 



W efffi) 



-fx 



(35) 



3(1 — fi m ,o) 

It is strongly supressed by small / 1; so this contribution can be sizeable only if the 
potential U{$>) is very steep today. However, steep potential would lead to the rapid 
acceleration of the scalar field, so the phantom phase would be very short in the past. 
Furthermore, the fast evolution of the scalar field together with large U\ = dU/d&(z = 0) 
would imply rapid change in time of the effective dark energy density. To elaborate on 
the latter point, let us consider the parameter 

dw eff 



Wi 



dz 



0). 



Observationally, \wi \ is not large: the WMAP analysis [T] gives —0.31 < w% < 1.12. On 
the other hand, making use of Eq. ^ and neglecting the terms with i*\ we obtain 



Wi = w hA + w hB , 



(36) 



W 1:A 



Wl, B 



+ 



Wi ; b are: 




1 




3(1 — f2 m) o, 




1 




3(1 — f2 m ,o, 


2 


1 




3(1 — Q m fi] 





(F3-6WO, 

[(1 + 5fi m , ) (F 2 - 6Q Ufi - 3Q mfi + 6) - 9Q mfi ] + 

9 



12F 2 + 



^ + 9F 2 



u,o 



301] [/n 



+24(1 - Q mfi ) + 6 



1 



'— F 2 + QVL Ufi + 3Q mfi - 6)' 



(37) 



(38) 



3(1 - n m , ) 2 

The parameter F 2 cannot be very large, see below, so for large u\ = Ui/<b' Q the value 



of w\ is controlled by w\^a term in the expression (36). To have sizeable contribution 



(35) and at the same time satisfy the observational constraint on w±, one would need 
the cancellation between F 3 and 6ui, which in turn would require strong fine-tuning. 
Barring this possibility, we arrive at the conclusion that 1, so the contribution 



(35) is very small. From now on we neglect it. 
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The remaining terms in Eq. (34) can be simplified by using (33): 

(i + / 2 )K) 2 



l + w. 



efffi 



(39) 



3(1 — tt m ,o) 

Thus, the phantom behaviour today is controlled entirely by j% = d 2 F/d§ 2 (z = 0). To 
have w e ff < — 1 at the present time, one requires that 

/ 2 < -1. 



Together with the bound (31), this implies that today the field $ must be close to a 



relatively sharp maximum of the function F(<&). Clearly, such a special state requires 
fine-tuning of both the function F(Q) and initial conditions in the theory. 



a. 



AF«1 



A$ «1 



FIG. 1. The shape of F($>). The present value of <3? must be near the maximum. 



We continue the discussion of the shape of F(Q) and recall the constraints on the 
time-dependence of the gravity constant. Given the small value of F\ and large value of 



Wbd,o, the constraint (16) is not hard to satisfy; note that this is in contrast to Ref. [20J. 



Much less trivial is the fact that the "cosmological" gravity constant has not changed 
much since BBN. It is clear from Q that G cosm {<&) is simply equal to F _1 ($). Barring 



strong cancellations, the constraint ( JlTj ) together with Eq. (10) imply that the variation 
of has been small since BBN, 



AF < 0.1 



(40) 
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and that F' is small at large z. Thus, the function F(<&) must have the shape shown in 

Fig. [g 

Let us now estimate the range of redshifts in which the dynamics of F is non-trivial, 
and the phantom effective equation of state can be realized. We do this by requiring 
that the value of F does not change much during this period. Since Fx is small, the 
second term in the expansion of F in redshift is relevant, and the estimate for maximum 
redshift is found from 

- FoZ 2 < AF 

2 2 max ~ s ' 



where the bound on AF is given in Eq. (40). Making use of Eq. (29) we get 

l(\M&o) 2 zL x < AF. (41) 
Let us denote by e the deviation of w e // from —1 today: 

1 + w e ff,o = -e- 



Using (39) we get the estimate (1 + /^X^o) 2 ~ 3e, and hence from (41) we find 

< 2AF\1 ± M 
3e |/ 2 | 

For reasonably strong phantom behavior (i.e., not very small e), ^ max is fairly small; 
roughly speaking, z max < 1. Note that similar result has been obtained within the 
reconstruction approach [7J [131 HTJ. At larger redshfts, F is frozen out, and the scalar 
field <3> reduces to quintessence. 

VI. NUMERICAL EXAMPLE 

Let us give a concrete example of a model with effective phantom behavior today 
and in the recent past. We note that the relatively large value of the parameter e = 
— (1 +w e ff i o) is obtained for fairly large $' , otherwise the maximum of F must be very 



sharp (i.e., / 2 must be large), see Eq. (39). So, we take, somewhat arbitrarily, $' = 0.3. 



We choose f2 mj o = 0.25, in rough agreement with observations, then Eq. (10) with Fi <C 1 



gives f2[/o = 0.73. Fairly strong phantom behavior, e = 0.1, is obtained with f 2 = —1.75. 



The estimate (42) then gives z max ~ 0.2. To satisfy all these requirements, we choose 



F(<&) as shown in Fig. [2j 
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0.99B 
0.996 
0.994 
0.992 
0.99 
0.9BB 

0.02 0.04 0.06 O.DB 0.1 <P 

FIG. 2. F($) in the numerical example. 

The shape of the potential U ($) is not constrained particularly strongly; our choice 
is shown in Fig. [3j 

2.2B 
2.27 
2.26 
2.25 
2.24 
2.23 

0.02 0,04 0.06 0.0B 0.1 <J> 

FIG. 3. U(&) in the numerical example 

With this choice, the effective equation of state depends on redshift as shown in 
Fig. |4j As expected, w e ff rapidly tends to —1 as z increases, and the field $ becomes 
indistiguishable from quintessence at z > 0.2. The field value does not change much: 
the change from redshift 0.2 to the present epoch is about A$ m 0.05. 

We have constructed a number of other examples satisfying the constraints of Section 
III; all of these examples have similar properties. 
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W'eff 

-1 
-1.02 
-1,04 
-1.06 
-1.0B 
-1.1 
-1,12 

0.15 0,3 

FIG. 4. w e ff as function of redshift. 

As we already pointed out, effective phantom dark energy requires both the special 
form of the function F(Q>) and fine-tuning of the initial value of the field $. To see 
the latter property explicitly, let us take the same functions and U($) as before 

and consider the evolution from redshift z = 0.75 to z = for different initial values of 
the field. If we vary the initial condition for $ within 15% interval around the central 
value yielding Fig. [4] (without varying the initial velocity $, for the sake of arument), 
the evolution changes considerably. In particular, the effective equation of state is as 
shown in Fig. [5j 



-1.1 
-1.2 

-13 
-1.4 
-1.5 

0,75 

FIG. 5. w e ff for initial conditions for $ deviating by ±15% from the central value. 
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Such a behaviour is not unexpected. The right choice of the initial value of the 
field ensures that the phantom phase begins in just right time, at some rather small 
redshift. For different initial values, the onset of the phantom behavior occurs at "wrong" 
redshifts, so one either has too large deviation of w e //,o from —1, or no deviation at all. 



VII. CONCLUSIONS 



In this paper we revisited the question of the possibility of the present phantom 
phase in scalar-tensor gravity. We have seen that it is possible to obtain and control 
effective phantom behavior even in simple scalar-tensor models, but this requires a lot 
of fine-tuning. First, the large present value of the Brans-Dicke parameter is obtained 
only if the scalar field $ is presently near the extremum of the function F determining 
the gravity constant. This is consistent with observable phantom property only if this 
extermum is a sharp maximum. Second, the small variation of the gravity constant 
since BBN requires that F flattens out at fairly low redshift. Finally, the whole picture 
is consistent with observations only for fine-tuned initial data. 

We conclude by noting that some of the unpleasant properties discussed in this paper 
must be present in scalar-tensor models for effective dark energy irrespectively of whether 
it is phantom or not. This remark applies, in particular, to the contrived shape of the 
function F and fine-tuning of initial conditions. Indeed, the fact that dF/dQ must be 
small today does not rely on the assumption of the phantom behavior. Furthermore, 
most of the analysis in Section V goes through provided that \d 2 F/d<& 2 \ is large enough 
at the present epoch, while the case \d 2 F/d<& 2 \ <C 1 corresponds to quintessence rather 
than genuine scalar-tensor gravity. All this makes scalar-tensor theory rather unlikely 
candidate for explaining the accelerated expansion of the Universe. 
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